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I discuss soft-gluon rcsummation and power corrections for event shape distributions, mostly in e + e~ annihilation. I 
consider specifically the thrust, the C parameter, and the class of angularities, and show how factorization techniques 
and dressed gluon exponentiation lead to predictive models of power corrections that are firmly grounded in pertur- 

■ bative QCD. The scaling rule for the shape function for angularities is derived as an example. Finally, I make a few 
i remarks on possible generalizations to hadron collisions, and on their relevance to LHC studies. 

o : 
<n : 

1. INTRODUCTION 

P . 
i-^ . 

\Q ' Quantum Chromodynamics is made special among phenomcnologically relevant field theories by the property of 
confinement. Given that the lagrangian is written in terms of fields that do not appear in the construction of the 
i true asymptotic states of the theory, it may seem surprising that perturbative calculations performed around the 
trivial vacuum have any relevance at all. The predictive power of perturbative QCD, in the presence of a kinematic 
scale Q 2 much larger than the confinement scale A 2 , is rescued by asymptotic freedom, combined with quantum- 
mechanical incoherence and gauge invariance. These are the necessary ingredients entering the proof of factorization 
theorems 1], which are the cornerstones of all PQCD calculations. 

Factorization theorems in essence provide a bound for the parametric size of nonperturbative corrections to high 
energy inclusive cross sections. Such corrections are typically suppressed by powers of the hard scale Q 2 . It should 

■ be emphasized, however, that factorization theorems are proven perturbatively, by examining the all-order struc- 
Oh, ture of long-distance singularities in Feynman diagrams. Their phenomenological relevance must then rely upon the 
O i additional (if plausible) assumption that confinement be a relatively soft process, happening without a violent rear- 
rangement of momentum configurations, as colored particles evolve away from the hard scattering event. Decades of 
experience with QCD phenomenology have taught us that this assumption is very well borne out by the data. Indeed, 

. j-h , granted this assumption, the bound on the size of nonperturbative effects provided by the factorization theorem is 
the first and simplest case of nonperturbative information extracted form QCD by purely perturbative methods. 

The general idea that perturbation theory, whenever genuine all-order information is available, can provide impor- 
tant clues to understand nonperturbative effects, has subsequently been applied successfully in a variety of situations. 
Typically, the perturbative expansion is found to diverge, and the uncertainty in the physical prediction due to this 
divergence is interpreted as a measure of the size of the expected nonperturbative correction. Specifically, the non- 
perturbative contribution must be ambiguous by an amount matching the uncertainty in the perturbative prediction. 
The assumption that the actual size of the nonperturbative corrections should be well represented by this ambiguity 
is sometimes referred to as ultraviolet dominance of power corrections Q ■ 

The all-order perturbative information required to begin any study of power corrections has mostly been provided by 
two complementary sources: renormalon-type calculations (reviewed in Q), which roughly speaking target running- 
coupling effects by summing up fermion bubble corrections to single gluon emission, and soft gluon resummations, 
which make use factorization and universality to compute leading multigluon contributions in the soft and collinear 
limits (for a recent review, see Q). Recently, it was shown that the two appoaches can be combined U yielding a 
strongly constrained and rather elegant model of power corrections in the Sudakov region. 

Event shape distributions in hard collisions are an especially interesting class of observables for power correction 
studies, and indeed a lot work has been done in the past several years on the subject, especially in the context of e + e~ 
annihilation and DIS |6|>Li- Event shape distributions, in fact, provide a continuous interpolation between processes 
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featuring mostly hard, perturbative radiation and configurations dominated by soft and collinear gluon emission. 
The corresponding theoretical prediction must then be constructed matching a variety of tools: NLO perturbative 
results for hard emissions, soft gluon resummations when the value of the event shape forces radiation to be soft, 
and finally models of power corrections very close to threshold. 

In general, models of power corrections involve nonperturbative parameters or functions, which must be determined 
from experiment, much as one does with parton distributions. The predictive power of these models must then rely 
on a degree of universality of soft radiation, which is well understood in perturbation theory, and must be assumend 
to hold nonperturbatively as well. By comparing theoretical predictions for different but related event shapes one 
can then test our understanding of QCD at or beyond the strict limits of applicability of perturbation theory. 

The application of these techniques has lead in recent years to quantitatively testable and quite successful models 
of power corrections. Here I will mostly discuss results obtained by Dressed Gluon Exponentiation (DGE) as 
applied to thrust @, |J, the C-parameter 0, and the class of angularities 0, 0. These examples show that 
current tools lead to simple, analytical, quantitative results that can readily be compared with experimental data. 
Most strikingly, leading power corrections to angularity distributions obey a simple scaling rule as a function of a 
continuous parameter, which gives a powerful test of our understanding of soft QCD in electron-positron annihilation. 
In Sect. © I will briefly summarize the formalism of shape functions for event shape distributions, and show how 
DGE provides a renormalon model for shape functions incorporating the constraints of NLL soft gluon resummation. 
I will use mostly thrust as a working example, comparing at the end with similar results obtained for the C parameter. 
In Sect. Js3 I wm discuss the class of angularities and derive the scaling rule, and finally in Sect. Q I will briefly 
comment on possible extensions of these techniques to hadron-hadron collisions. 



2. SOFT GLUON EFFECTS FOR EVENT SHAPE DISTRIBUTIONS 

An event shape distribution is a weighted cross section, assigning a prescribed value to a specific infrared and 
collinear safe combination of the momenta of final state particles in a high energy collision. In the case of e + e~ 
annihilation, let F m (pi 7 . . . ,p m ) be one such combination, computed for an m-parton final state. The distribution of 
the associated event shape / is then 

% = ^ E / dLIPS ™ l-^™! 2 6 (/ - F m(Pl, ■ ■ ■ ,Pm)) , (1) 

where M m is the appropriate matrix clement. In the following, I will consider event shapes / which vanish in the 
limit of a pencil-like two-jet event. A prime and well-known example is r = 1 — T, with T the thrust, 



T = max 



Ei Wi 



(2) 



which I will use below to illustrate the general features of the approach. Other shapes I will consider include the 
C-parameter, 

C-3 3 V (Pi-Pi? (3) 
which does not require a maximization procedure, and the one-parameter class of angularities, 

where rapidity rji and transverse momentum px,» are computed with respect to the thrust axis, so that for a = one 
verifies that r = r. 

The common feature of these event shapes, which opens the way to an all-order perturbative analysis and to studies 
of power corrections, is the fact that for small values of / all radiation is constrained to be soft or collinear. As a 
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consequence, the distributions develop double logarithmic singularities of Sudakov type, which can (and must) be 
resummed thanks to the universal properties of soft radiation and to the factorizability of the cross section in the 
Sudakov limit. In QCD, resummation displays the ambiguity of perturbation theory, originating from the presence 
of the running coupling evaluated at soft scales. This leads naturally to models of power corrections. I will now 
illustrate the general features of the method using thrust as an example. 



2.1. Resummation 

In order to resum singular contributions to the thrust distribution in the limit r — > 0, one needs to take a 
Laplace transform, which factorizes the (5-function constraint fixing the value of r. Logarithmic contributions then 
exponentiate according to 



aire 



. 1 da 
a dr 



exp 



/'^(e— -l)(B(a s (uQ 2 ))+2 H %A(a s { q 2 )) 

Jo u \ Juzqz <r 



u 2 Q 2 1 



(5) 



The pattern of exponentiation is highly nontrivial, since the double logarithms of the ordinary perturbative expansion 
turn into single logarithms in the exponent. Generically one finds a structure of the form 



iry pc 
13] 



2 k 



Ck P L p — > exp Lgi(a s L) + g 2 (a s L) + a s g 3 (a s L) 



(G) 



where L is the logarithm of the transformed variable, L = log v in this case. Leading logarithms (LL) are generated 
to all orders by the function g\, which is completely determined by the knowledge of the anomalous dimension 
A(a s ) to one loop. Next-to- leading logarithms (NLL), corresponding to the function g 2 , require the knowledge of 
A{a s ) to two loops and B(a s ) to one loop. NLL accuracy is the common standard for resummation of event shape 
distributions. Note however that Sudakov resummation, expressed here by Eq. (J5J, in general involves one more 
function, D (a s {u 2 Q 2 )^ . This function is associated with wide-angle soft gluon emission, and is process-dependent, 
unlike the anomalous dimension A. To any finite logarithmic accuracy, the contributions of D can be reproduced by 
modifying B in a process-dependent manner, and for the event shapes discussed here D does not give any contribution 
at NLL level. The fact that D dpends on the scale uQ, however, has important consequences on power corrections, 
as discussed below. 



2.2. Power Corrections 

One can deduce from Eq. J^J) (where the integration variable u in the exponent plays the role of r) that for small 
values of r there are two relevant momentum scales: tQ 2 and t 2 Q 2 . This can be understood from the physical 
picture underlying Sudakov factorization: at small r gluon radiation can be organized into jets of particles collinear 
to the primary partons, with invariant mass proportional to \/tQ, plus the contribution of wide-angle soft gluons, 
characterized by their total energy tQ. It is natural to expect that power corrections will be organized by these 
two scales, and thus be of the form (A 2 /(rQ 2 )) m and (A 2 /(r 2 Q 2 ))™ respectively. Clearly, when r ~ A/Q all power 
corrections of this second kind become important, and must be collectively taken into account. Power corrections in 
the larger scale, on the other hand, become important only when r ~ A 2 /Q 2 , a value which is too small to be relevant 
for LEP fits. The need for power corrections is apparent in Eq. (JSJ), since the integrals over momentum scales are 
perturbatively ill-defined because of the Landau singularity of the running coupling. An elegant way to summarize 
the nonperturbative information encoded in Eq. (JSJ was described in 0] . The basic assumptions are the applicability 
of the factorization underlying Eq. all the way down to values of r such that A 2 ~ t 2 Q 2 <C tQ 2 <C Q 2 , and the 
existence of a nonperturbative definition of the running coupling rendering the scale integrals well defined. Consider 
then, for example, the term in Eq. (JSJ containing the anomalous dimension A. In order to disentangle perturbative 
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and nonperturbative domains, one can simply introduce a factorization scale /i, switch the order of the q and u 
integrations, and define 



I" 1 fin f uQ do 2 

S^Q 2 ) ee / _(e- u " -1) / A r A{a s (a 2 )) = S NP (v/Q,») S P t(", Q, m) , 

*,<„/<,,„) . £' f A (a ,(, 2)) £° | - 1) - f i (i)" A„,, 2 ) + O (f) . (7, 



The last equality expresses a set of nonperturbative contributions to the Sudakov exponent in terms of moments of 
the anomalous dimension A at low scales. These moments, 



A„(M 2 ) = i T dq 2 q n - 2 A(a s (q 2 )) , 
71 Ja 



(8) 



are not computable in perturbation theory: much like parton distributions, they should be measured for a given 
observable at a given factorization scale, and then used to predict different observables, based on their universality 
properties. In general, the full set of leading nonperturbative corrections will involve also moments of the function D, 
which also parametrize power corrections of the form (Ais/Q) n . These corrections thus have a universal component, 
expressed in terms of the anomalous dimension A, and a process-dependent component given by the function D. At 
power accuracy, it is natural to disentangle the contributions of B and D by requiring that the function B should 
be the same appearing in the resummation formula for DIS structure functions, where the corresponding D function 
is known to vanish [15j. Expressions like Eq. Q provide a framework to test universality, or to construct specific 
models of power corrections. To summarize the effects of the parameters A n (/i 2 ) one can use them to build up a 
"shape function" , according to 



exp 



SnpOVQ.m)! = / &e-^/«/ T)NP (e,/*). (9) 



Here e can be interpreted as the total energy carried into the final state by soft gluons at scales below /i. Confining 
oneself to the leading power correction, corresponding to the first moment Ai(/x 2 ), one recovers the result of the 
"tube model" 16]: that nonperturbative effects shift the distribution away from the small r region by an amount 
proportional to the average energy carried away by soft radiation. Subleading moments provide additional smearing. 



2.3. Dressed gluon exponentiation 

The shape function idea is very general, and can be used both to test universality, by connecting power corrections 
to related event shapes 01 > or to construct models based on factorization and Lorentz invariance in specific cases 0] . 
One can get more detailed predictions by making stronger assumptions: for example, one can apply a renormalon 
model, and study the corresponding power corrections in the Sudakov region. This is the basic idea underlying 
dressed gluon exponentiation (DGE) 5]. I will now summarize the basic steps of this method, using thrust as an 
example. 

First of all, one computes the single gluon contribution to the event shape under study, for a gluon of nonvanishing 
virtuality £ = k 2 /Q 2 . This is the characteristic function of the dispersive approach [l9j, |20( to power corrections. 
Since we are interested in the Sudakov region, we need to retain only terms that are singular as r — > 0. Given the 
characteristic function r), one can write a clean representation of the single gluon cross section by introducing 
a Borel representation for the strong coupling. One defines 

A{£,Q 2 ) = r du ^ u (Q 2 /^T U — c KU ■ ( 10 ) 
Jo ' wu 

This amounts to an analytic continuation of the strong coupling at the scale k 2 from the euclidean to the timelike 
region, formally valid in the large-n/ limit. The factor e KU is renormalization-scheme dependent, with n = 5/3 in the 
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MS scheme. The single dressed gluon cross section is then 



1 da fn . rfi\ Cf 
o~ dr 2p 







s 1 nOO 

A{tiQ 2 ) = %L du (Q 2 /A 2 y u B(t,u), 
z Po Jo 



(11) 



where in the last equality I introduced the Borel function B(t,u), which contains the physical information on the 
thrust distribution. The strategy of performing all integrals except the one over the Borel parameter yields at the 
end a trasparent representation for power corrections. In the case of thrust, the terms responsible for logarithmic 
enhancements in r) are given by 



log 



which gives a Borel function of the form 



B(t,u)1 = 2c h 



sin nil 



-l-2u -1-u 
- T — T 



1 — u 2 — u 



(12) 



(13) 



Note that already at this stage one can make several useful observations. Poles in B(t,u) at, say, u — uo, would 
correspond to renormalon singularities in the distribution, and expected power corrections of size (A/Q) 2u °; in fact, 
B(t, u) has no such poles: the would-be singularities at u = 1,2 are cancelled by the factor sin7ru, and poles at u = 
cancel between the two terms in square brackets, because of the infrared safety of thrust. Renormalons arise when 
taking moments of the distribution, because of the convergence constraints on the Borel integral at small values of 
t. One also notes that the first term in Eq. (|13l) is associated with soft wide-angle radiation, since it generates in 
Eq. i|ll[l terms proportional to (A/ (tQ)) 2u . Similarly, the second term in Eq. (|13fl is associated with the jet function, 
contaning collinear as well as soft enhancements. 

The key step in DGE is to note that at LL level Sudakov resummation yields a simple exponentiation of the one- 
gluon emission cross section in moment space. One can then retain all large-n/ information, and the corresponding 
model of power corrections, in the Sudakov exponent by simply using the single dressed gluon cross section as kernel 
of exponentiation. One defines 

,fc+ioo dv 



1 da\ 

where the Sudakov exponent is now given by 



fc-i 



27ri 



E{v, Q 2 ) 



/ dr(l 
Jo 



1 da\ 



e vr exp Q 2 )] , 



B T (v,u) . 



(14) 



(15) 



, a dr /SDG 

Here the single dressed gluon cross section is defined by Eq. l]l II) , virtual corrections have been taken into account by 
subtracting the value of the Laplace transform at v = 0, and in the second equality the Borel function B T (v, u) for the 
Sudakov exponent has been defined. As usual, all integrals are performed except the one on the Borel parameter u. 
Although formally similar to Eq. (jl 1|> for the single dressed gluon cross section, Eq. I|15J1 has a much richer physical 
content, displayed by the nontrivial renormalon structure of the Borel function B T {v 1 u). This is a consequence of the 
fact that exponentiation, subject to the constraint of energy conservation, has promoted the single gluon result to a 
genuine approximation for multigluon emission 1 . In the specific case of the thrust, the result for the Borel function 
is 



B T (v,u) = 2e K 



sin ttu 



T(-2u) {v 2u - 1) - - r(-«) K - 1) - + 



2 -it 



(16) 



1 In fact, the exponent E(u,Q 2 ) has a natural interpr etat ion in terms of the Borel representation of Sudakov anomalous dimensions 
associated respectively with soft and collinear radiation Il5l . Corrections subleading in change the the logarithmic behavior of the 
cross section as well as the size of the residues of poles in the Borel plane, but they are not expected to modify the analytic structure of 
B T (y,u), which determines which power corrections are actually present. 
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It is useful to compare this result with the one for the C-parameter [lOj 




r(-«) („«-i)^ + 



1 



1 



u 



+ 



2 -it 



1 



) 



(17) 



Clearly, Eqs. and l|17(l are very similar; there are, however, important differences, which highlight the degree of 
universality to be expected in comparing the two distributions, and which can in principle be tested experimentally. 
First of all, the Borel functions contain perturbative information on Sudakov logarithms: although the exponentiation 
was performed assuming independent gluon radiation, it can be shown Q that one can upgrade the formalism to 
NLL accuracy by simply replacing the running coupling Eq. I|10|) with the two-loop expression, and by changing 
renormalization scheme, including in the constant k the contribution of terms singular as x —> 1 in the NLO Altarclli- 
Parisi splitting function (the "gluon bremsstrahlung" scheme 21]). Beyond NLL, the coefficients of all subleading 
logarithms can be computed in the large n/ limit, by simply expanding the Borel function in powers of it, and 
replacing u n — > n!(6oa s /7r)™ +1 . Computing subleading logarithms uncovers the factorial growth of their coefficients, 
and can be used to gauge the reliability of perturbative resummation in different kinematical regimes. Next, one 
may observe that the infrared safety of r and C is once again reflected in the cancellation of the poles of both Borel 
functions at u = 0. One also notes that wide-angle soft radiation and collinear gluons contribute as before two 
separate terms, and the jet functions (the terms proportional to F(— u)) are identical for the two observables 2 . They 
contribute renormalons at u — 1, 2, corresponding to exponentiated power corrections of the form (A 2 v/Q 2 ) 1 ' 2 . Soft 
gluon contributions, on the other hand, have renormalons at u = m/2, for all odd values of to, yielding the leading 
power corrections (Av/Q) m , and they are quantitatively different for r and C, distinguishing the two observables. 
Specifically, by taking the ratio of the two "soft functions" (the terms proportional to T(—2u)) and expanding in 
powers of u, one can verify that the two observables begin to differ perturbatively at NNLL level, as predicted 
in [22j | . but the growth of the coefficients of further subleading logarithms is weaker for the C-parameter than for 
the thrust. Similarly, if one boldly takes the large-n/ residues of the poles of the Borel functions as a reasonable 
estimate of the size of the corresponding power corrections, one observes that (Ai//Q) m corrections are systematically 
smaller for the C parameter: the two shape functions should therefore differ, and one expects that the resummed 
perturbative prediction, as well as the approximation of the shape function by a constant shift, should work better 
phcnomenologically for C than for r. 



The discussion in Sect. @ illustrates the predictive power of DGE. Another interesting application concerns 
angularities, whose definition, Eq. can be rewritten for massless particles as 



with 0i the angle with respect to the thrust axis. Angularities (so christened in |l2j) were introduced in |23( as 
auxiliary shape variables used to tame nonglobal logarithms |2^| for observables related to out-of-jet energy flow. 
They have several remarkable features, which make them very interesting for our understanding of QCD at the edge 
of the perturbative domain . First of all, they are characterized by a tunable parameter a, which can be used to 
interpolate between different shapes, or to bring to focus specific momentum configurations in a continuous way. The 
parameter a must satisfy a < 2 for infrared safety, and a tighter restriction a < 1 is required in order to preserve 
a relatively simple resummation in the Sudakov region, in the form of Eq. (jSJ): for 1 < a < 2 further logarithmic 
singularities associated with jet recoil must be taken into account [2£j. For a = 1, one recognizes that t\ — B, the 



2 Note that B c (u,u) is computed for a rescaled C parameter, c = C/6 



3. THE CLASS OF ANGULARITIES 




(18) 
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broadening; for a = 0, tq = 1 — T; for negative a, events dominated by high rapidity give increasingly suppressed 
contributions to r a , which in turn suppresses power corrections of collinear origin; finally, for a — > -co the distribution 
becomes a 5 function at r a = 0, with a strenght given by the total cross section. 

Remarkably, although the relative weights of rapidity and transverse momentum change with a, it is possible to 
derive a resummation formula [2^ of the form of Eq. (|SJ), valid for a < 1. Indeed, at NLL accuracy one can write 



1 r ^2 
hx[a(u,a)] = B(a s (uQZ))(e-^ 2/(2 -" } -l)+2 £ ^ ^ A (a s (q 2 )) (e-^/Qr _ 1) 



(19) 



The a dependence of Sudakov logarithms is clearly nontrivial: as an example, the function g\{a s L) responsible for 
leading logarithms in Eq. JSJl is given by 



gi(x,a) 



4 2 - a 



/3q 1 — a x 



1- x 

2- a 



In (1 - x) 



2 - a 



In 1 



2 - a 



(20) 



Power corrections, however, turn out to have a much simpler a dependence |ll| . Performing the analysis leading to 
Eq. JSJ), one easily finds that in the nonperturbative region all moments of the anomalous dimension A are multiplied 
by a simple common factor 



a! q V) 



1 1 



1 — a n 



dq\^A(a s {q 2 )) = T^-AgV) 
1 — a 



(21) 



As a consequence, the Laplace transform of the shape function defined by Eq. © obeys a simple and remarkable 
scaling rule 



fa 



NP 



fa 



NP 



Q 



, K 



l/(l-a) 



(22) 



which should be experimentally testable without great effort using existing LEP data. Lac king a direct experimental 
analysis, the scaling rule was tested against the output of PYTHIA, with positive results [ll|,|2i|. The physical picture 
underlying the scaling rule is appealing, and once again reminiscent of the "tube model" . The relevant feature of the 
radiation pattern is boost invariance, which applies to soft gluons emitted in the two-jet limit, since such emissions 
are correctly represented by the eikonal approximation. This means that soft gluons contribute to the event shape a 
rapidity- independent amount, and in turn the integration over rapidity simply measures the size of the region where 
gluons can be emitted without strongly affecting the event shape. This region scales with (1 — a) -1 . The derivation 
of the scaling rule neglects correlations between gluons emitted into the opposite hemispheres defined by the thrust 
axis, which however are expected to become important only for a > 1, the region excluded by the present treatment. 
The main assumption is then that nonperturbative soft radiation should share the property of boost inveriance with 
the relatively harder perturbative component which is treated by resummation. This is the nonperturbative property 
that would be directly tested by an experimental study of Eq. I|22|l . 

One can go further, and study subleading power corrections, mostly related to radiation collinear to the hard jets. 
Applying the same method to the anomalous dimension B, and to the subleading terms generated by A, once again 
one finds a simple pattern: the Laplace transform of the cross section can be expressed in factorized form, introducing 
a subleading shape function, as 



a (y, a) = a PT (z/, k, a) / aiN p 



Q' 



9a, NP 



Q 2 



-, K 



(23) 



Clearly, as a becomes large and negative, collinear power corrections are expected to become more and more negligible, 
and the scaling rule Eq. (|22|) is expected to hold with increasing accuracy. 

In the light of the discussion of Sect. (0), it is interesting to verify whether these nice features of angularities are 
preserved in specific models for the shape function, such as DGE. This is particularly relevant in this case, since the 
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key result, Eq. (|22|l . is related to boost invariance, which is broken by the formal introduction of gluon virtuality, 
which is a necessary tool of renormalon analysis. The study of angularities with DGE was performed in [l^. It 
is technically nontrivial, since the interplay of the parameter a with gluon virtuality and phase space constraints 
makes it difficult to extract the a dependence analytically. The first step is to find an appropriate generalization 
of the definition of angularity to the case of single massive gluon emission. Such definition should have the correct 
limit as £ — > 0, reduce to known results for thrust as a — > 0, and be simple enough to keep the computational task 
manageable. At one loop, one such definition is 



(1 



\l-a/2 



(1 _ x _ £)l-a/2 (l _ Xk+ f) -/a + {] „ k) 



(24) 



where x n — 2p n -Q/Q 2 are the customary energy fraction variables, and the definition applies to the phase space region 
where the gluon is soft; xt is then the (anti)quark energy fraction. With this definition, it is possible to construct 
the Borel function for the Sudakov exponent for angularities, in analogy with Eqs. (|f 6[> and (|f 7fl . Remarkably, the 
soft component of the Borel function, responsible for all leading power corrections, is just the expected rescaling of 
Eq. (J3J, 



1 



„„, sin7ru , , , o, \ 2 
2 e r(-2u) (v 2u - 1) - 

ITU V ' U 



(25) 



which leads once again to the scaling rule, Eq. Ij22(l . Collinear power corrections are much more difficult to handle 
analitycally, and the collinear counterpart of Eq. (|25|) can at best be expressed in terms of a one-dimensional integral 
representation, which reduces to combinations of hypergeometric functions for rational values of a. This is however 
enough to classify the singularities in the Borel parameter u, and thus the pattern of power corrections. One indeed 
finds that all these subleading power corrections can be organized in a single shape function g, depending only on 
the combination v/Q 2 ~ a , as in Eq. (|23|) . DGE thus confirms the general scaling behavior found from resummation, 
showing that the introduction of gluon virtuality does not spoil the effects of boost invariance in the Sudakov limit. 
Thrust, jet masses, angularities and the C parameter are all found to have closely related pattern of power corrections, 
highlighted by the scaling rule relating generic angularities to the thrust. Clearly, this is a highly predictive framework, 
and our understanding of soft radiation in the two-jet limit can be put to stringent tests. 



4. HADRON COLLIDER EVENT SHAPES 

As we approach the expected date for the start up of the Large Hadron Collider at CERN, it is natural and 
appropriate to ask whether tools like those described here could be applicable in the environment of hadron collisions, 
and, if so, to what extent power corrections might be relevant to our understanding of the data, at the extreme energies 
available at the LHC. Beginning with the second question, one might naively observe that A/Q must be a very small 
number for any reasonable value of the hard scale Q that one might envisage at the LHC. ft would not be wise, 
however, to neglect power correction studies on this basis, for at least three different reasons. 

First of all, it has already been shown at the Tevatron that power corrections have an impact even on observables 
which are largely dominated by high pj_ events, and even at very high energy [2t| . In Ref. [27| . Mangano considered 
the single-jet inclusive E± distribution, comparing data at different CM energies. One sees that ratios of cross sections 
at different energies do not have the proper scaling behavior dictated by NLO QCD, but the correct behavior can 
be recovered by including a power correction determined by a single parameter associated with the normalization of 
the jet transverse energy. The reason is that even a small shift in the jet E±_ is amplified in the distribution by the 
fact that the cross section is falling steeply for increasing E±, so that 5a /a ~ —nSE± if a ~ One can expect 

that in general power corrections will be important for the determination of jet energy scales, and in turn accurate 
knowledge of these scales may well turn out to be crucial for many high energy studies. 

The second reason is that almost any LHC observable will require, before it can be compared with a theoretical 
prediction, a subtraction of all hadronic activity unrelated to the hard scattering, loosely referred to as "underlying 
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event". There is currently very little, if any theoretical control on the underlying event, which at the LHC will 
contain a mix of multiple parton scattering, beam-beam interactions and soft radiation associated with the selected 
hard process. Even Monte-Carlo methods have difficulties in finding the proper tuning to describe this kind of 
physics (see, for example, |28|). In this context, the lesson of power correction studies in the gentler environment of 
e + e~ annihilation or DIS is that we might learn to discriminate between the different components of soft radiation in 
hadron collisions. On the one hand there are soft and collinear gluons associated with the hard scattering event: their 
effects arc in principle computable in PQCD, using generalizations of the known techniques, and their distribution in 
phase space and in the space of color configurations will be nontrivial and predictable. On the other hand, there is 
soft radiation which is in practice out of reach for the techniques of PQCD, such as minijets due to multiple parton 
scattering or soft gluons arising from beam-beam interactions. This second kind of radiation fills phase space with 
a high degree of uniformity, and will have to be modelled with different techniques, including Monte-Carlo tools. 
This kind of statistical modelling is bound to be more successful if we can first achieve a better understanding of the 
"pure" hard scattering process, including the energy and color flow that it generates at all scales. 

Finally, it should be emphasized that the physics of event shapes at hadron colliders is interesting for its own 
sake, as a probe to understand hadronization in terms of both momentum and color flow. It is well known that for 
hadron collisions, where most processes involve four partons already at Born level, Sudakov resummation is expressed 
in terms of anomalous dimension matrices that tie together color exchange and momentum flow. This interplay is 
bound to influence the pattern of power corrections as well, and studies of this kind may well lead to deep insights 
into the mechanics of color neutralization and hadron formation. 

Preliminary studies of rcsummed event shapes at hadron colliders have already been performed [29! l3fj| , and a 
generalization of the notion of angularity to a hadronic environment has been proposed |3l| . For most of these 
proposals, a primary concern is that of suppressing the contributions of particles close to the beam axis, where 
beam remnants interfere with all measurements. It should be noted, however, that soft radiation associated with 
the underlying event tends to fill phase space, including regions separated in rapidity from both the beam and the 
high-p^ jets. It would therefore be of great interest to find event shapes designed to focus on this kind of wide-angle 
radiation, where it would be most useful to disentangle soft gluons generated by the hard scattering from the genuine 
underlying event (a pioneering study with a similar goal is [32l| ) . A promising avenue of investigation might be the 
use of observables such as those introduced in [2^], joint distributions correlating energy flow in a chosen angular 
region £1 with a standard event shape such as ordinary angularity. The form of such a correlation is 



a e 2 , a) = - J2 \ M ( N )\ 2 ~ MN)) S(e 2 - r^(N) - r^(N)) , (26) 



where fn(N) is an observable related to energy flow into the angular region f2, away from hard jets, while Ta (N) 
are the contributions to angularity from the two hemispheres defined by the thrust axis. 

In a hadronic environment, one might envisage measuring angularities with respect to the current jet axis, as 
suggested in |3lj . or even introducing a third parameter £3 to constrain radiation near the beam remnants. Tuning 
the various parameters and energy threshold in such observables one would be able to focus on different regions of 
phase space, forcing the observable to be more or less inclusive for soft gluons, or for particles collinear either to 
the beam or to the current jets. Notice that boost invariance, underlying the scaling rule for angularities, is lost in 
hadron collisions, where correlations between beam jets and current jets must be taken into account. 

Clearly, event shape studies at hadron colliders are in their early days. I believe that such studies will be both 
instrumental to further our understanding of QCD, and very helpful in order to exploit the full potential of the LHC, 
a task which will require a solid understanding of strong interactions as much as good skills in the building and 
testing of new physics models. 
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